Generic model of an atom laser 



OO 

< 

O 

(N 



I 

o 
o 



> 
00 

o 

00 



X3 
O 

o 



X 
J3 



B. Kneer,i'ld T. Wong,^ K. Vogel,i W. P. Schleich,i and D. F. Walls^ 
^ Abteilung fiir Quantenphysik, Universitdt Ulm, D-89069 Ulm, Germany 
^Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand 
(Submitted to Physical Review A: June 18, 1998; revised manuscript submitted: August 20, 1998) 

We present a generic model of an atom laser by including a pump and loss term in the Gross- 
Pitaevskii equation. We show that there exists a threshold for the pump above which the mean 
matter field assumes a non-vanishing value in steady-state. We study the transient regime of this 
atom laser and find oscillations around the stationary solution even in the presence of a loss term. 
These oscillations are damped away when we introduce a position dependent loss term. For this 
case we present a modified Thomas-Fermi solution that takes into account the pump and loss. Our 
generic model of an atom laser is analogous to the semi-classical theory of the laser. 
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I. INTRODUCTION 



With the recent experiments on Bose-Einstein conden- 
sation [p]-p^ an atom laser, that is a device which pro- 
duces an intense coherent beam of atoms by a stimulated 
process p^ , p^ , has become feasible. Already, the MIT 
group has realized a pulsed atom laser and has pro- 
vided evidence for the process of coherent matter-wave 
amplification in the formation of a Bose condensate [0 . 

How to describe the matter field of an atom laser? Rate 
equations are simple. But they cannot answer this ques- 
tion since they do not contain any coherence. In contrast 
a microscopic and fully quantum mechanical treatment 
can answer the question but is not easy to handle. What 
we need is a theory that includes coherence but is still 
simple. In the present paper we develop such a theory — a 
generic model of the atom laser. 

The optical laser has three essential ingredients: (i) 
a resonator for the electromagnetic field, (ii) an atomic 
medium, and (iii) an excitation mechanism for the atoms. 
We start the laser cycle by preparing the atoms in an 
excited state using this excitation mechanism, which in 
general is incoherent. The radiation emitted by the atom 
amplifies the electromagnetic field in a mode of the res- 
onator. To be efficient the frequency of the mode has 
to match appropriately the frequency of the transition. 
The boundary conditions set by the resonator determine 
via the Helmholtz equation the spatial part of the mode 
function. The ultimate goal is to transfer the excitation 
of the atom into a macroscopic excitation of the field 
mode. In this way we transfer the energy used to excite 
the atom via the gain medium into coherent excitation of 
the field mode. In order to make use of the radiation we 
have to couple it out of the resonator. We compensate 
for this loss by continuously re-pumping the medium. 

The goal of an atom laser is completely analogous: We 
want to create a macroscopic coherent excitation of a 
mode in a resonator for atoms. Hence, in an atom laser 
the atoms play the role of the field excitation of the op- 
tical laser. Since atoms cannot be created or annihilated 



the means to achieve lasing are different. Indeed, in an 
atom laser there is no "real" laser medium: It is the 
same atom that goes through the laser process. We only 
manipulate the internal degrees of freedom and center- 
of-mass motion of the atoms. In particular we want to 
force their center-of-mass motion into a specific quantum 
state of the resonator. The resonator for the atoms is a 
binding potential such as provided by a trap. The spa- 
tial part of the mode function of this atomic resonator 
follows from the time independent Schrodinger equation 
]l8[ . Moreover, we focus on the ground state of the trap. 
As in the optical laser we want a macroscopic excitation 
of this mode, that is we strive to have as many atoms as 
possible in one quantum state. This is the phenomenon 
of Bose-Einstein condensation. As in the optical laser we 
need to couple the atomic wave out of the atomic res- 
onator. In order to have a continuous wave (cw) atom 
laser, we have to continuously feed in more atoms. 

There exist two different approaches towards a theo- 
retical description of a cw atom laser: The first one relies 
on rate equations [p^-pT| whereas the second one derives 
a quantum mechanical master equation ||2^-p7[|. In the 
present paper we suggest a third approach which makes 
heavily use of the close analogy between an atom laser 
and an optical laser. In the latter case it turned out that 



a classical treatment of the electromagnetic field [£8 29 
was sufficient to describe many features of the laser. Can 
we therefore devise a semi-classical theory of the atom 
laser? 

The semi-classical laser theory replaces the electro- 
magnetic field operator E{r,t) for the field inside the 
laser cavity by the expectation value £{r,t) = {E{r,t)). 
The equation of motion for £ is the wave equation of 
Maxwell's electrodynamics driven by the polarization of 
the laser medium. An additional term introduced phe- 
nomenologically takes into account the loss of the cav- 
ity. The polarization of the laser medium follows from a 
microscopic, quantum mechanical description of the in- 
ternal structure of the atoms. Quantum mechanics rules 
the atoms whereas classical Maxwell's wave theory deter- 
mines the electromagnetic field. These are the essential 
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ideas of semi-classical laser theory. 

The semi-classical laser equations do not prefer any 
particular phase. Nevertheless by choosing an arbitrary 
phase we can describe many properties of the electromag- 
netic field. Furthermore, we have to start with a non- 
vanishing seed field in order to obtain a non-vanishing 
solution for the electromagnetic field with a fixed phase. 

In our model of an atom laser we replace the matter- 
wave field represented by a field operator ^(r, i) by a 
scalar mean field = (^'(r, i)). This is analo- 

gous to replacing the field operator E by its expectation 
value £ in semi-classical laser theory. The well-known 
Gross-Pitaevskii equation [^,0 plays now the role of 
Maxwell's wave equation. It defines the equation of mo- 
tion for tp{r,t). Similar to the semi-classical theory of 
the optical laser we have to break the symmetry of the 
equation of motion for ^ in order to have a non- vanishing 
value for -ip. For massive particles this is more problem- 
atic than for photons since they cannot be created or 
annihilated, and for all quantum states with fixed parti- 
cle number we have ('5) — 0. Nevertheless the concept 
of spontaneously broken symmetry turned out to be very 
useful to describe properties of a condensate, in partic- 
ular interference effects. For more detailed discussions 
see Refs. In contrast to the driven electromag- 

netic wave equation the Gross-Pitaevskii equation does 
not contain a gain term analogous to the polarization. 
This reflects the fact that in Bose-Einstein condensation 
there is no "medium" in the trap. We therefore add a 
phenomenological pump term. Moreover, as in the elec- 
tromagnetic case we have to add a loss term. 

Despite the similarity there is a fundamental differ- 
ence in the two equations of motion. In the absence of 
a medium Maxwell's wave equation is linear. In con- 
trast the Gross-Pitaevskii equation is non-linear. This is 
a manifestation of the interaction of atoms. 

The crucial part of any laser is the stimulated ampli- 
fication process. Different mechanisms for matter-wave 
amplification have been suggested and discussed: optical 
cooling |l9[-[23|, elastic collisions by evaporative cooling 
p5| , p6|j33| 7dissociation of molecules ||2^ , and cooling by a 
thermal reservoir p^ ]. For different schemes of pumping 
a condensate we refer to Refs. ||35[-^. Our model does 
not rely on a specific mechanism, but can be adapted to 
any mechanism where Bose enhancement is present. 

The paper is organized as follows: In Sec. ||, we gen- 
eralize the Gross-Pitaevskii equation by including gain 
and loss terms. One important loss is due to coupling 
the atom wave out of the resonator. However, similar 
to semi-classical laser theory where usually the field in- 
side the laser cavity is investigated, we restrict ourselves 
to the matter-wave field inside the resonator. We there- 
fore do not go into details of an output coupler [^s]-^. 
We present the stationary solutions of our equations and 
perform a stability analysis. We find a threshold behav- 
ior similar to the optical laser. Moreover, we calculate 
the time dependent solution and show how it converges 



to a quasi-stationary solution. In general, the quasi- 
stationary solution we find does not coincide with the 
time independent solution but shows some oscillatory be- 
havior around it. These oscillations disappear when we 
modify our equations in Sec. Ill by introducing a space 
dependent loss. This loss can be thought of as a conse- 
quence of collisions between condensed and un-condensed 
atoms at the edges of the condensate. For this improved 
model we present a modified Thomas-Fermi approxima- 
tion for the stationary solution. Section IV summarizes 
our results. 



II. ELEMENTARY MODEL 

In this section we summarize our generic model of an 
atom laser. 



A. Formulation of the model 

In particular, we add the phenomenological gain and 
loss terms to the Gross-Pitaevskii equation (GPE). In 
this way we couple the GPE to an equation governing 
the number of un-condensed atoms. 

The Gross-Pitaevskii equation (GPE) for the mean 
field ip — ip{r) of the condensed atoms of mass m in 
a trap potential V{r) reads HEsl 



at 2to 



(1) 



The non-linear term C/qIV'PV' with C/q = Airfi^as/in takes 
into account two-particle interactions where denotes 
the s-wave scattering length. 

The GPE has been very successful in describing the 
properties of Bose-Einstein condensates. However, in the 
present form it cannot describe the growth of or the loss 
of atoms out of a condensate. Indeed, the GPE keeps the 
number of atoms 



Nr. 



(2) 



constant. 

In order to overcome this problem we introduce two ad- 
ditional terms in the GPE: a loss term and a gain term. 
The loss term 



(3) 



leads to an exponential decay of the number of atoms in 
the condensate. In contrast the gain term 



ih 



(4) 



leads to an increase of atoms in the condensate. Here 
is the number of atoms outside the condensate, that is, 
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the un-condensed atoms, and T is the rate for the transi- 
tion of these atoms into the condensate. We regard this 
as a generic pump mechanism of an atom laser since 
it contains ah'eady the Bose enhancement as we will see 
later. 

When we add the loss term and the gain term (^) 
to the GPE (|l|), we arrive at the generalized GPE 



ot 2m 



(5) 



For the number of un-condensed atoms we assume 
the rate equation 



(6) 



The first term reflects the fact that we generate with a 
rate i?„ atoms in an un-condensed state from an infinite 
reservoir of atoms. The term —^uNu takes into account 
that atoms can escape from our system without being 
trapped in the condensed state. The last term describes 
the transition of the atoms into the condensate. Since it 
is proportional to the number of atoms Nc in the con- 
densate, it contains the Bose-enhancement factor. 

In this way we have coupled the generalized GPE gov- 
erning the condensate to the rate equation governing the 
number of un-condensed atoms. These two equations (|^) 
and (^) are the foundations of our model. 

The number Nc of condensed atoms follows from the 

fieralized GPE (||) with the help of the definition, Eq. 
. From Eq. (||) we obtain the rate equation 



Nc = TNuNc - -fcNc 



(7) 



for the number of atoms in the condensate. 

We conclude this section by noting that rate equations 
similar to Eqs. ( [6[) a nd (Q) have already been discussed 
in the literature |19| , pO| . However, in the present paper 
we replace the rate equation by the generalized GPE 
(^. This equation obviously contains the rate equation 
but in addition the coherence. 



B. Stationary solutions and stability analysis 

Before we turn to the stationary solution of the gener- 
alized GPE, Eq. §), we first discuss the stationary solu- 
tions of the rate equations (^) and (|^). 



1. Rate equations 



The rate equations 
tionary solutions: 
The solution 



and (|^) have two possible sta- 



N^""^ = 0, 



(8) 



is reminiscent of the optical laser below threshold where 
the intensity of the laser vanishes. 
The other stationary solution 



Mis) _^_lu 



(9) 



corresponds to the laser above threshold. 

We now perform a stability analysis of these solutions. 
For this purpose we introduce small deviations 

nu{t) = Nu{t)-Ni'\ 

nc{t) = Nc[t)~Ni'\ (10) 

and arrive at the linearized equations 

fin = -(7n + r7Vi^))n„ - TNi'\c. 

nc = -{lc-TNi'^)nc + TN'^c'\^. (11) 

A stability analysis of these equations shows that the 
stationary solution in Eqs. (p|) is stable for 



Ru < R 



th _ 7c7« 



(12) 



Likewise, our stability analysis shows that Eqs. (^ are a 
stable stationary solution for 



Ru > 



Iclu 



R 



th 



(13) 



Therefore, there exists a threshold i?*''. When the pump 
rate i?„ is below the threshold, the number N,, of atoms 
in the condensate vanishes. When it is above, it is non- 
vanishing. In that case the steady-state number of atoms 
in the condensate grows linearly with the pump rate Ru- 



2. Generalized GPE 

We now turn to the discussion of the stationary solu- 
tions of our generic model of the atom laser, that is of the 
generalized GPE coupled to the rate equation for Nu- 

The stable stationary solution of Eq. (^ corresponding 
to a vanishing number of atoms in the condensate is 



V'^"^ = 0. 



(14) 



This is the stationary solution of the mean field below 
threshold. 

Above threshold the stationary solution ■0^*' of the 
generalized GPE (||) is identical to the stationary solution 
of the conventional GPE. Indeed, when we substitute the 
ansatz ^ = exp(— into Eq. (||) and note that 
the gain and loss terms cancel each other as a result of 
the stationary solution, Eq. (0), we arrive at 
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^^(«) = - A^M + y(r)?A(^' + C/o|?/'("'|V^''- (15) 
2m 

Here, /x denotes the chemical potential. 

A stability analysis of this equation leads to collec- 
tive excitations which in general cannot be treated 
analytically. For such an analysis in the case of a 
one-dimensional harmonic oscillator we refer to the Ap- 
pendix. 



C. Transient behavior 

In the preceding section we have discussed the steady- 
state solutions of both the rate equations (^) and (0) and 
the matter-field equations (^) and of the atom laser. 
In the present section we address the question if and how 
the matter field evolves into a stationary state from an 
initial condition. Here we focus on the case above thresh- 
old. Since the equations are non-linear we solve them 
numerically for the case of an one-dimensional harmonic 
oscillator potential 



V{x) = —muj^x^. 



(16) 



This one-dimensional trap potential already shows the 
essential features of the atom laser. 

We therefore analyze the one-dimensional generalized 
GPE 



., 9 , , , 



2m dx^ 
i , I 



1 



-hVN,, 



(17) 



describing the matter field, coupled to the equation 



(18) 



for the number of un-condensed atoms. We use the split- 
operator method to solve the generalized GPE. 
This technique has already been used successfully for the 
ordinary GPE [Q. Moreover, we supplement this nu- 
merical analysis by an analytical treatment of the solu- 
tion in the long-time limit where we apply the method of 
Rcf. to one dimension. For the details of the analyt- 
ical approach we refer to the Appendix. 



1. Initial conditions and parameters 

The semi-classical theory of the optical laser cannot 
explain the initial start-up of the laser since it does not 
contain spontaneous emission: A non- vanishing seed elec- 
tromagnetic field starts the laser. Similarly, in our model 
of the atom laser we start from an initial condition for 
the mean matter field iplx^t — 0) that is non-zero. This 
is the seed field for our atom laser. 



Let us illustrate this by considering for the moment 
the "natural" initial condition ^{x,t = 0) = and 
iV„(0) 0. This implies iVc(O) = 0. We therefore start 
without any atoms in the system. Moreover, we con- 
sider a pump rate above threshold. This choice of 
initial conditions leads to the unstable solution Nc = 
and Nu = Ru/lu- Indeed, any small perturbation in 
leads to a completely different behavior and Nc and Nu 
approach the stable solutions N^f'^ — Ru/lc — 7ti/r and 



Nu^'^ = 7c/r, as discussed in Sec. H B . Therefore, we use 
the different initial condition Nc{Q) ^ 1 for our numeri- 
cal simulations [^. We keep the condition A^ti(O) = 0. 

For our numerical calculations we take the parameter 
Ux/{tujoa) = 0.008, where a = ^JhjimLj) is the width of 
the ground state of the harmonic oscillator. As in con- 
ventional laser theory p8| , p9[ , where one usually has a 
high-Q cavity, we require 7c w. Since in the steady 
state we do not want to have too many atoms in the 
non-condensate part, we require 7c ^ P, as suggested by 
Eq. (^. In order to have a reasonably small threshold 
value we choose 7„ appropriately, as indicated by Eq. 
( p^ ) . We also need to choose Ru such that we have a suf- 
ficiently large number of atoms in the condensate, that 



Ruhc - 7«/r » 1. 



2. Average properties of the matter field 

In Fig. 0(a) we show the time dependent solutions 
Nc{t) and iV„(i) of the rate equations (^) and (|^). 
For very small times the atoms accumulate in the un- 
condensed phase: Due to the small number of atoms in 
the condensate the Bose enhancement is not yet effective 
and they cannot make a transition into the condensate. 
As soon as we have a significant number of atoms in the 
condensate, rapidly approaches its stationary value 



7c/r. Additional atoms then essentially end up 



in the condensate where the number of atoms slowly ap- 
proaches its stationary value ivi"'' = Ru/"fc — 7m/P- 

For the numerical analysis of the Eqs. ( p^ ) and ( p^ ) 
defining our generic model of the atom laser, we use the 
same parameters. For the initial condition 'ip{x,t — 0) 
of the generalized GPE we use the quantum mechanical 
ground state of the trap, normalized in such a way that 
A^c(O) = J\4>{x,t — 0)\^ dx. In order to get some feel- 
ing for the accuracy of our split-operator technique for 
the generalized GPE we calculate Nc{t) = J \il;{x,t)\'^ dx 
which according to Sec. HA has to coincide with the so- 



lution of the rate equations. 

We gain insight into the time dependence of |i/'(a;,i)P 
by calculating its first moment 



it) = 



and second moment 



/ x|V'(a;,t)p dx 

J \ij{x,t)\'^ dx 



(19) 
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J\ij(x,t)\'^dx ■ 



(20) 



We note that the symmetry of the trap and the gen- 
eraUzed GPE ensure that the symmetry of our initial 
condition is preserved, that is 



it) = 0. 



(21) 
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FIG. 1. Generic model of an atom laser: transient regime, 
(a) The number Nc of atoms in the condensate for the 
spatially independent (solid line) and spatially dependent 
(dash-dotted line) loss, and the number of un-condensed 
atoms Nu for the spatially independent (dashed line) and spa- 
tially dependent (dotted line) loss approach their steady-state 
values 10^ and 1, respectively. The scale on the vertical axis is 
logarithmic. The inset shows the early stage of the transient 
regime. Here, we cannot distinguish the curves corresponding 
to the different spatial damping, since for small atom numbers 
the spatial dependence does not play a big role, (b) The width 
Ax /a of the mean field as a function of time displays collec- 
tive excitations with period T = 3.57/uj as apparent from 
the inset. For the spatially homogeneous damping the oscil- 
lations are un-damped whereas for the spatially dependent 
loss they are damped away as shown in (c). The parameters 



10" 



10" 



are Ru/uj = 10^ F/lu = 10"^ 7„/lj 
and Ux / (Jiuocl) = 0.008. As initial conditions we have used 
iV„(0) = 0, iVc(O) = 10"^ and for ij){x,t = 0) we have used 
the quantum mechanical ground state of the trap, normalized 
in such a way that J \ijj{x,t = dx = Nc{Q)- For (c) we 
have used a spatially dependent loss 7c (a;) described by Eq. 
(E) with -y'c/uJ = 0.192, xo/a = 5, and a/a = 1. 



Hence the width 



Ax(t) 



-x^ = ^Jx^t) 



(22) 



of |V'(a;,t)P is governed by the second moment, only. 

In Fig. |^(b) we plot the width Ax. We note that Ax 
increases as a function of time and approaches a steady- 
state value for large times. Moreover, it oscillates around 
this steady-state value. In the inset we magnify these os- 
cillations. We note that the oscillations do not decay 
Therefore Ax does not approach a time indepen- 
dent value and there is no steady state in a strict sense. 
We hence refer to this solution as the quasi-stationary 
solution. 

In the Appendix we derive an analytical expression 
for the frequency ^2 of these oscillations and find ^2 = 
V3w ^ 1.73w. From Fig. |(b) we read off f^a = 1.76w 
which is in good agreement with our prediction. 



3. Matter field 

The oscillations in the width of the distribution are a 
manifestation of collective excitations of the conden- 
sate. Indeed, the whole distribution oscillates as shown 
in Fig. Ij. Here we display the mean field for a time inter- 
val where the numbers of atoms in the condensed state 
and the un-condensed state have already reached their 
steady-state value. 

In order to study this oscillatory behavior in more de- 
tail, we compare in Fig. |^(a) the solution of the time 
independent GPE ( p^ in one dimension to the quasi- 
stationary solution of Eqs. (^^ and (^8|). Here we have 
depicted |V'(a;,t)p for time moments where the width is 
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FIG. 2. Quasi-stationary solution |'!/'(a;,f)p of our 
atom-laser equations for a time interval when the number of 
atoms in the condensate has already reached its steady-state 
value. We have used the same parameters as in Fig. nl 



maximal, average, and minimal. Indeed, the solution of 
Eqs. (^7|) and (jlj) oscillates around the solution of the 
time independent GPE [|5|-|7)- 




III. IMPROVED MODEL 



The collective excitations discussed in the preceding 
section depend on the initial condition ijj{x, t = 0). This 
is due to the fact that so far there is no mode selection 
mechanism in our model which would favor one station- 
ary solution of the GPE over the others. This fact is 
similar to the optical multi-mode laser [^,^. We can 
accomplish a single-mode atom laser when we allow for a 
space dependent loss. This is the topic of the present sec- 
tion where we formulate an improved model for an atom 
laser. 



A. Formulation of the problem 

Real losses in a trapped condensate are spatially de- 
pendent: For example, a possible loss mechanism are 
collisions with un-condensed atoms which are preferably 
located at the edge of the condensate, see Refs. 
Further, atoms can be lost at the edge of the conden- 
sate because of a finite trapping potential. Moreover, 
the goal of a matter-wave output-coupler is to create a 
directed coherent atomic beam. All of these reasons re- 
quire a spatially dependent loss term. In contrast to this 
loss term we keep the pump term spatially independent. 
This reflects the fact that we assume that the condensate 
is pumped from a cloud of cold thermal atoms which is 
larger in size. 

We therefore consider a modified generalized GPE 



-5 x/a 5 

FIG. 3. Quasi-stationary and stationary solutions of our 
atom-laser equations. On the top (a) we compare the solu- 
tion of the time independent Gross-Pitaevskii equation (dia- 
monds) to the quasi-stationary solution of Eq. (hTI) for times 
where the width is maximal (dashed line) , average (solid line) , 
and minimal (dash-dotted line), respectively. These times 
correspond to Ljt = 1000.54, 1001.45, 1002.35. On the bottom 
(b) we consider the case of spatially dependent loss. We depict 
(not to scale) the loss rate ^c{x) by a dotted line. In this case 
collective excitations are damped out and in steady-state the 
solutions of our atom laser equations ( p^ ) and ( |25[ ) coincide 
with the time independent ground-state solution (diamonds) 
of Eqs. ( |28[ ) and (^^. In both figures we have used the pa- 
rameters of Fig. Py. 



where the decay rate 7c(r) is spatially dependent. Hence, 
the rate equation for the number Nc of condensed 
atoms now translates into 



TV, = TN,,Nr 



(24) 



at Zm 



(23) 



Note that due to the space dependent loss rate we cannot 
express the integral in Eq. ( p4| ) by Nc- Therefore we do 
not have a rate equation for Nc anymore. However, the 
rate equation 
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iV„ = Ru - luNu - TN.Nu 



(25) 



2. Lasing threshold 



for the number of un-condensed atoms remains the 
same. 

Equations ( p3| ) and ( p5[ ) are the two fundamental equa- 
tions of our improved model of an atom laser. 



B. Stationary solution 

We now study the stationary solution of Eqs. (|2^) and 
(^). We first derive the steady-state expressions and 
then discuss the laser threshold in the presence of a po- 
sition dependent loss term. We conclude by introducing 
a modified Thomas-Fermi approximation. 



1. Matter field in steady-state 



When we follow an analysis similar to the one of Sec. 



II B , we find two sets of stationary solutions of the Eqs. 
(If), (H), and (H). The first set with 



Ru 

lu 

0, 




(26) 



corresponds to the solution below threshold. 
The second set 



U 



7c(r)|^J^)(r)|2(iV 



Ru 



lu + TN, 



Ru 



lu 



/7,(r)|^(^)(r)|2d3^ 



(27) 
(28) 



corresponds to the solutions above threshold as we show 
now. Here, ~ \Ncf^\~^l'^'\\)^'^^{r) denotes the sta- 

tionary solution of Eq. ( P3| ) normalized to unity. This 
solution is defined by the time- independent modified gen- 
eralized GPE 



= _ AV'(^) + V (r) + J/olV'^'^lV^'^ 
2m 

- Ift7,(r)v;(^) + ^ftr7v„^(^\ 



(29) 



following from Eq. (|23| ) using the ansatz iiij^t) — 
exp(— i/ii/fi,)-0'^''^(r) where ^ denotes the "chemical po- 
tential." 

fs) Is) 

These expressions for Nu and Nc are not explicit 
since they in volve t he stationary solution ■0^*' (r) of Eq. 
(p9|). In Sec. tllBS we derive an approximate analytical 
expression and in Sec. IIIB5 we find a fully numerical 
solution for -i/;^''-' (i")- 



A stability analysis of the above solutions amounts to 
calculating the collective excitations of the modified gen- 
eralized GPE d^). This is only possible numerically. We 
therefore pursue a strategy where we solve numerically 
for the full time dependence of ?/'(r, t). In order to gain 
some insight into the threshold condition we first discuss 
simple physical arguments. 

These considerations rely on the fact that the number 
of atoms cannot be negative. With the help of Eq. ( [2^ ) 
we then determine the laser threshold 



R 



th 



1 



f /7c(r)|v4^^(r)|^rf 



2 



by setting Njf^ = 0. Close to threshold the number Nc 
of atoms in the condensate is close to zero. We can there- 
fore neglect the non-linear contribution in the modified 
generalized GPE and arrive at a linear Schrodinger equa- 
tion with pump and loss. With a position dependence of 
the loss term that favors the normalized ground-state en- 
ergy solution 4>Q (r) of the linear Schrodinger equation, we 

can replace ip'^^ (r) by 0o (r) in Eq. (pO|) and arrive at 



(30) 



As) 



R'^ = ^ j 7,(r)|</>o(r)|2d3r. 



(31) 



When the loss shape is flat around the localization of 
0o(r), that is, around the center r = of the trap, we 
can factor out 7c (0) from the integral in Eq. ( pl| ) and 
find the lasing threshold 



R 



th 



luiM 



(32) 



We conclude this section by noting that in the limit of 
a space independent loss t he abo ve results reduce to the 
corresponding ones of Sec. II B I. 



3. Modified Thomas-Fermi solution 

In this section we derive an approximate but analytical 
expression for the stationary state of the modified gen- 
eralized GPE. In the case of the conventional GPE it is 
the so-called Thomas-Fermi (TF) approximation which 
describes the steady state of the condensate [Q. The 
phase of the stationary solution of the GPE as well as 
the phase of the TF approximation is constant. How- 
ever, for a position dependent loss it turns out that the 
phase of the stationary matter-wave field also depends 
on the position. A spatially dependent phase leads to a 
non-vanishing current of the condensate. We expect this 
feature to be important for the coherence properties of 
a cw atom laser. Therefore, we introduce a modified TF 
solution. 

We start from the time independent form of the modi- 
fied generalized GPE, Eq. (09). For the present problem 
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a hydrodynamic treatment is more convenient. We there- 
fore consider the density 



and velocity 



p(r)E^|^(«)(r)|^ 



v(r) = {h/Tn)V(t>{r) 



(33) 



(34) 



of the condensate. Here 0(r) is the phase of the mean 
field following from the ansatz 



^Z^^'') (r) = VpWe"^^''^- (35) 
Indeed, with Eq. (p9|) we find the equations 



r- 

2m 



-/i + -rnv^ + y (r) + U^p 



and 



2 = [-V • V - 7, (r) 



(36) 



(37) 



We make the assumption that the density profile is slowly 
varying. Since we are interested in loss shapes where the 
loss at the center of the condensate is uniform, these as- 
sumptions are reasonable at the center of the condensate. 
However, at the edges of the condensate where the change 
of the loss is large, we expect this assumption to break 
down. This is the idea of our modified Thomas-Fermi 
approximation. 

This assumption allows us to neglect the terms V^y^ 
and Vy^ in the two equations. However, in order to ful- 
fill the second equation we have to retain the derivative 
V • V in velocity. In this approximation we arrive at the 
approximate expression 



p(v) 



1 



/i — -mv^ (r) — V (r) 



for the density and the differential equation 



V • V (r) = -7c (r) 



(38) 



(39) 



for the velocity. 

Since the density cannot be negative, these expressions 
are only valid for that volume V of space where 



p(r) cx /X - (r) - V" (r) > 0. 



The shape of this volume V depends on the potential 
V{y) and the loss rate 7c (r) via the velocity v(r). 

Note, that the velocity v and consequently the phase 
(p depends on the rate F, the number of un-condensed 
atoms, and the shape of the loss rate. Indeed, it is only 
the spatial dependence of the loss rate that determines 
the spatial profile of the velocity v(r). The "kinetic en- 
ergy" of the condensate plays a role similar to the trap 
potential in shaping the density. 



The expressions, Eqs. ( |38D and (^9|), for the density 
and velocity are not explicit. Indeed, the chemical po- 
tential p, is still a free parameter. Moreover, the number 
of un-condensed atoms Nu'' is coupled to Eqs. ( ^ and 
(H) and is given by Eq. (|^). 

We therefore have to solve the various constraints in 



a self-consistent way: 
atoms reads 



The number Nr of condensed 



N^/) = / p{v) d^r. 



(41) 



On the other hand this quantity N, 
and reads 



is) 



is given by Eq. (|2S 



/v7c(r)pi(r) d^r 



(42) 



is) 1 

where pi{v) — [Nc ]^ p(v) is the density normalized to 
unity. Note, that V is defined by the condition that the 
density is non-negative, Eq. (^0|). 

The number of un-condensed atoms then follows from 
Eq. ( ^7|) . Hence we have to solve the three Eqs. (40), (|4l]), 
and (^) for the three unknowns: the modified Thomas- 
Fermi volume V, the chemical potential /i, and the num- 
ber Nr of atoms in the condensate. 



4- Usual Thomas-Fermi solution 



The most important consequence of our position de- 
pendent loss term is a non-vanishing velocity and there- 
fore a spatially dependent phase. How does this compare 
to the usual Thomas-Fermi approximation? Here we con- 
sider the density 



p{v) = 



M - ^(r) 



(43) 



(40) and the velocity 



V = 0. 



(44) 



We first note that this ansatz is in contradiction to Eq. 
(^). Nevertheless we can use it to investigate how the 
kinetic energy potential influences the density p. Similar 
to the last section we insert our ansatz into Eqs. ( ^ ) to 
( ^ ) and find the Thomas-Fermi volume V, the chemical 

potential p., and the number of atoms Nc ■ 
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5. Exact solution 

In this section we adapt the numerical methods de- 
veloped to find the ground-state solution of the ordi- 
nary GPE to the present problem. This method 
evolves the wave function for a fixed atom number in 
imaginary time. The solution is normalized to unity after 
each time step. The evolution in imaginary time atten- 
uates the differences between the arbitrary initial wave 
function and the ground-state solution of the ordinary 
GPE. 

In our model of an atom laser, Eq. (^), we have gen- 
eralized the GPE by a pump and a loss term. Therefore 
the number of atoms is not fixed but is governed by Eq. 
(p5|). In order to find the ground-state solution of our 
time independent Eqs. (p8| ) and (|2^), we evolve an ar- 
bitrary initial wave function according to our modified 
generalized GPE, Eq. ( |23|), in imaginary time, using the 
split-operator technique [ [48p9| ] . After each time step we 
normalize the wave function to unity. Then we update 
the number of atoms with the help of Eq. (^8|) and use 
it for the next time step. We repeat this procedure until 
the wave function and the atom number has converged 
to a stationary value. This method finds the stationary 
state of our improved model of an atom laser in a self- 
consistent way. The most important result is that the 
pump and spatial loss gives a space dependent phase to 
the stationary mean-field ip'-^^r). 



C. Results 

As in Sec. || we now specify the potential V{r) to be 
the one-dimensional harmonic oscillator potential, Eq. 
(|l6|). Moreover, we study two models of a position de- 
pendent loss rate: a sum 

of two Gaussians and a sum 



7c (x) = 7^ 



{x + Xo) + CT^ {x - xo) +0-2 



(46) 



of two Lorentzians. 

In both cases ct is a measure of the width of the cor- 
responding distributions and 7^ is the maximal loss rate 
when the two curves have negligible overlap. In the next 
section we choose the locations ±xo of the maximal loss 
such that they sit where the density of the condensate 
falls off. 



due to the presence of the space dependence in the loss 
. For this analysis we use the Gaussian loss term, Eq. 
(p5|). We adjust 7^ such that we obtain the same final 
number of atoms in the condensate as in the spatially 
independent case. The width of the loss is comparable 
to a typical scale, such as the width of the ground state 
of the harmonic oscillator. 

In Fig. |^(a) we compare the time evolution of the num- 
bers of condensed and un-condensed atoms for the posi- 
tion dependent loss with those of the position indepen- 
dent loss. We see that the overall behavior is not that 
different. Therefore our original rate equations of the el- 
ementary model still approximate very well the time evo- 
lution of Nc and following from the improved model. 

We then show in Fig. 0(c) the scaled width of the mean 
field. Indeed, the position dependent loss damps out the 
collective excitations, and a steady state is reached. We 
emphasize that this is true for any initial non-vanishing 
mean field. 

In Fig. ^(b) we compare and contrast the stationary 
mean field to the numerical solution of the time inde- 
pendent equation. They coincide with each other. In- 
deed, for many initial conditions and set of parameters we 
have noticed that the mean field obtained by evolving the 
modified generalized GPE over a sufficiently long enough 
time such that the transient oscillations are damped out 
is equivalent to the time independent ground-state energy 
solution of the modified generalized GPE. This holds as 
long as the loss is confined to the edges of the condensate 
|]62| . This is the case in the present example as shown in 
Fig. ^(b) where we also display the loss function in ar- 
bitrary units. The overlap between the stationary mean 
field and the loss function is essential for the number of 
atoms in the condensate as is apparent from Eq. (p8|). 

When we compare Fig. ||(a) and (b), we note that the 
modulus of the steady-state mean field in the presence 
of a loss located at the edges of the condensate is very 
similar to the ground-state energy solution of the conven- 
tional time independent GPE, Eq. (p^. However, this is 
not true for t he ph ase of the mean field as already dis- 
cussed in Sec. IIIB. 



A spatially dependent loss term located at the edges of 
the condensate damps out the collective excitations. In 
our theory of an atom laser this plays the role of mode 
selection, i.e., the ground-state energy solution survives 
whereas the excited states, i.e., the collective excitations, 
are damped away due to the fact that their overall spread 
in position space is greater than that of the ground state. 



2. Stationary solution 



1. Transient behavior 

First we want to show that the collective excitations 
which emerged in the elementary model are damped out 



(a) Modified Thomas-Fermi solution. We now turn to 
the discussion of the modified Thomas- Fermi solution de- 



rived in Sec. HI B 3 for a general potential. In the present 
section we restrict this analysis to a one-dimensional har- 
monic oscillator potential, Eq. (p^). Moreover, we choose 
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FIG. 4. Comparison between the two Thomas- Fermi and 
the numerical solutions of the atom laser equations in the 
presence of space dependent loss. We display density (up- 
per part) and current (lower part) of the numerical (solid 
line), modified TF (dashed line), and usual TF (dash-dotted 
line) solution. On the top we also show the space depen- 
dence of the loss term 70(2;). Note the position of the peaks, 
the height is not to scale. Here we have used the parame- 
ters R-a/uj = 5 X 10^, V/lu = 0.5, 7„/cj = 1, 7c/lj = 5, and 
Ux/{huja) 5z 0.008. In the left column the loss is concentrated 
at a; = ±5a, whereas on the right it is at a; = ±3a. Moreover, 
the widths are different a/a = l (left) and a/a = 0.5 (right). 

the Lorentzian loss rate, Eq. (^). This choice is solely 
motivated by the fact that we can perform the resulting 
integral of the differential equation (|39| ) . 

We start by first summarizing the Eqs. (^) to (|4|) for 
the chemical potential ^, the Thomas-Fermi radius R, 



and the number of condensed atoms N, 



-muj^R^ = 0, 



R 



p{x) dx = 0, 



(47) 
(48) 



Ru 



!^r1c{x)pi{x) dx 



(49) 



determining the modified Thomas-Fermi solution for the 
one-dimensional harmonic trap and loss. From Eq. (|3^ ) 
we obtain the one-dimensional density 



n mv (x) muj X 



(50) 



We find the velocity v when we substitute the Lorentzian 
loss rate, Eq. (|4^), into Eq. and integrate which 

yields 



V (x) = TNuX — 7^(7 



arctan 



arctan 
- Xo 



Xq 



a 



(51) 



Here we have set w(0) = in order to preserve the 
symmetry of the solution of the mean field, Eq. ([5^), 
p[x) — p{—x). The latter holds because the harmonic 
oscillator potential, the spatially dependent loss, and the 
generalized GPE show this symmetry. 
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(h) Usual Thomas- Fermi solution. In order to get a 
feeling how in our one-dimensional model the density p 
is influenced by the velocity v we also discuss the usual 
TF solution. For the one-dimensional harmonic trap the 
density, Eq. reads 



p{x) = — 



-^2 2 

a — m—LU X 
^ 2 



and the velocity, Eq. Q), reads 

V{x) EE 0. 



(52) 



(53) 



When we use this usual TF approximation, we easily find 
from Eq. (^) with v = the TF radius R = y/2jijrnu^. 
We then use Eq. ( ^8|) to calculate the number of atoms 
in the condensate as a function of the chemical potential 
fi. Inverting this equation we can express the chemical 
potential p = [3Nc'^^Ux\^uj/{aV2)]^/^ as a function of 
the atom number A^i"-*. Finally, we use Eq. ( ^ ) to de- 

( s) 

termine the number Nc ' of condensed atoms, the only 
unknown quantity. We emphasize that only this last step 
has to be done numerically. 

(c) Discussion. Figure ^ shows the modified TF, the 
usual TF and the fully numerical solution. We have cho- 
sen two different shapes of the Lorentzian loss curve: In 
the right column the width of the individual Lorentzians 
is half the size of the one on the left column. Moreover, 
we note that their location is different. We depict the 
modified TF solution by dashed curves, the usual TF so- 
lution by dash-dotted ones, and the numerical solution by 
solid lines. Again we show in arbitrary units the shape 
of the loss 7c (x) by the dotted curve. 

The density is shown in the upper part of Fig. |^. We 
see that for the parameters of the left column both TF ap- 
proximations for the density work quite well. The mod- 
ified TF solution approximates the center better than 
the usual TF solution. However, the radius derived from 
the usual TF approximation is a better estimate of the 
edge of the condensate than the modified TF radius. For 
the set of parameters used in the right column of Fig. ^ 
both TF approximations for the density do not work that 
well anymore. The reason for this break-down is that 
the condensate reaches too far into the loss region where 
the density varies strongly and the derivatives neglected 
in the derivation of Eqs. (|3|) and ( p^ become impor- 
tant. Nevertheless our modified TF solution shows at 
least qualitatively the same behavior as the fully numer- 
ical solution. Surprisingly, at regions towards the peaks 
of the loss where the loss rate is at its highest, the density 
is first decreased and then increased in the modified TF 
over the usual TF solution. This is understandable from 
Eq. (^8|) when we consider in the lower part of Fig. |^ the 
current j{x) = p{x)v{x). We note that at these regions 
first the current which is proportional to the velocity is 
at its highest and then decreases. 

We see that for the current, i.e. the velocity, the mod- 
ified TF solution is a very good approximation around 




'0 R^/uj 10^ 

FIG. 5. Comparison between analytical, but approximate, 
and numerical, but exact solutions of the atom laser equa- 
tions. For this comparison we analyze the threshold behavior 
of the atom laser: Condensate population Nc"'^ in steady-state 
as a function of the pump strength _R„. Over a wide range 
of pump strengths both approximate solutions, the modified 
(solid line) and the usual (dashed line) TF solution accu- 
rately describe the fully numerical solution (diamonds). In 
the region close to threshold the two approximate solutions, 
which are almost identical, clearly deviate from the fully nu- 
merical solution as shown in the inset in the upper left cor- 
ner. The other inset magnifies a region far above thresh- 
old where the modified TF solution is a better approxima- 
tion than the usual TF solution. The parameters used are: 
T/lu = 10"^, 7„/a; = 10, j'^/lj = 1, a/a = 1 and xo/a = 10, 
and Ua:/{hiua) ^ 0.008. 

the central region. It gives a good estimate for the overall 
behavior of the current (velocity) for both sets of param- 
eters. Around the central region of the condensate the 
current is approximately a linear function of the position, 
illustrating the flow of atoms to the ends where they are 
predominately lost from the peaks. At the tails of the 
condensate the modified TF solution even predicts a ve- 
locity changing the direction. This is in contrast to the 
fully numerical solution. The usual TF solution has zero 
velocity by default. 

We solve numerically the Eqs. (^ ) to (^9|) for the mod- 
ified Thomas-Fermi solution to find the condensate pop- 
ulation for various pump strengths and display this as 
the solid curve in Fig. ^. The dashed curve corresponds 
to the usual TF solution. We display the fully numerical 
solution with the help of diamonds. 

Both approximate analytical curves agree quite well 
with the fully numerical solution over a wide region of 
different pump strengths. However, at pump rates just 
above threshold the agreement is not that good which is 
shown in the upper left inset of Fig. ^ . Both the usual 
TF and the modified TF curves lie on top of each other 
and over-estimate the numerical result. The fully numer- 
ical solution crosses the horizontal axis at = 201 
whereas the two TF solutions cross at = 198. Close 
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to threshold, the spatial extent of the two approximate 
wave functions is becoming smaller with decreasing atom 
number. They are predominately located around the rel- 
atively flat loss region and the losses appear to be position 
independent. Therefore the laser threshold for the two 
approximate solutions is predicted by Eq. (^2|). In con- 
trast, the shape of the numerically calculated mean field 
can vary, and is in fact of Gaussian form close to thresh- 
old. Hence our approximation of the threshold, Eq. (31), 
agrees rather remarkably with the numerically calculated 
threshold. The deviation of the approximate threshold, 
Eq. (jH) , reflects the fact that even close to threshold the 
loss function does not appear spatially constant to the 
mean field for the parameters of Fig. When we recall 
that the TP approximation is not good at low atom num- 
bers it is rather surprising to see such small differences 
between the TF solutions and the numerical solution. 

The lower right inset of Fig. |5| zooms in on a region of 
higher pump strengths far above threshold where the spa- 
tial structure of the loss plays a role. Here the modified 
TF solution approximates the fully numerical solution 
(diamonds) better than the usual TF solution. This is 
understandable since the modified TF solution takes into 
account pump and spatially dependent loss by allowing 
for a spatially dependent phase, i.e. velocity. 



IV. CONCLUSIONS 

In summary, we have constructed a theory of an atom 
laser that is analogous to semi-classical laser theory. The 
matter-wave equation is a generalized Gross-Pitaevskii 
equation with additional loss and gain terms. We de- 
rive the lasing threshold and describe the build-up of the 
coherent mean field of a condensate. 

The elementary model uses a spatially homogeneous 
loss. Here we find un-damped collective excitations. 
Therefore the final mean field depends on the initial mean 
field: The known stationary ground state of the GPE 
which is the desired lasing mode cannot be reached in 
general. 

The improved model has a natural mode selection built 
in by a space dependent loss. In this way we achieve the 
desired single lasing mode. 

We have derived a modified Thomas-Fermi solution 
for the steady-state mean field. This solution takes into 
account the effects of a pump term and a position de- 
pendent loss term. In contrast to a constant phase of 
the usual Thomas-Fermi solution, the modified Thomas- 
Fermi solution has a spatially dependent phase, i.e. ve- 
locity, due to the permanent fiow of atoms in and out of 
the condensate. The modified Thomas-Fermi solution is 
a good approximation in regions where the loss shape is 
slowly varying and for sufficiently large atom numbers. 

We emphasize that our model of an atom laser is very 
simple and rather general. Therefore, we can apply it to 
different experimental configurations of cw atom lasers 



|63| , or current experiments |0-|l|l provided the evapora- 
tive cooling process (boson amplification) and the loading 
of the trap is run continuously, and an output coupling 
mechanism is applied. 
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APPENDIX: COLLECTIVE EXCITATIONS IN 
ONE DIMENSION 

Collective excitations are usually discussed in three di- 
mensions [ ^Tj . Within the framework of the Thomas- 
Fermi approximation, Stringari |5l|| has calculated an- 
alytically the excitation spectrum of a condensate in a 
three-dimensional isotropic harmonic trap. Since our nu- 
merical solution of Eq. (^) is done for a one-dimensional 
harmonic trap, we apply Stringari's method to a one- 
dimensional harmonic trap of frequency lu. We substi- 
tute 



V'(x,t) = Vp(M)e**("'*^ 



(Al) 



into Eq. (|^) and obtain after some algebra the hydrody- 
namic equations 



and 



| + |:(H + (7e-riv.)p = o 



dv d /m 2 TT 



(A2) 



2m dx dx^ 

for the "density" 

p{x,t) = 'tj)*{x,t)'4){x,t) 

and the "velocity" 

h d(j){x,t) 



v{x, t) 



dx 



(A3) 



(A4) 



(A5) 



Note, that we have also introduced the chemical potential 
/i which is space independent. 

We now consider small deviations 
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Sv = V — u*-*^ — w, 

n„ = iV„ - Njf^ ^Nu- 7c/r (A6) 

of p, V, and iV„ from their stationary values w*-*-* = 0, 

and A^u"' = 7c/r- Furthermore, we make the Thomas- 
Fer mi ap proximation, that is, we neglect the last term in 
Eq. ( |A3| ) and approximate the stationary solution p^'^\x) 
by the Thomas-Fcrmi solution pQ] 



p^^>{x) 



(A7) 



where the chemical potential p is defined by the normal- 
ization integral. We then arrive at the linearized equa- 
tions 



dt 



and 



<5/9+^ (p(^)w) -Tp'^'lr 







(A8) 



(A9) 



We combine these two equations to eliminate v and arrive 
at 



Tp'^^'hu. 



(AlO) 



For large times, when Nu{t) and Nc{t) have already 
reached their stationary value, we can neglect the inho- 
mogeneous t erm T p^'^'^hu- 

Equation ( A10| ) cannot be solved without knowledge 
of the potential V{x). We restrict ourselves to the case 
of a harmonic trap, that is. 



V(x) 



2 2 

—muj X . 
2 



(All) 



In order to solve Eq. ( A10| ), wc introduce the scaled vari- 
able ^ — x/R, where R = \j2plmuP- is the Thomas- 
Fermi radius of the condensate. Using the ansatz 



6p{x,t) = Asm{flt + Lp)y{x/R) 



(A12) 



we obtain an ordinary differential equation for y(^) which 
reads 



d_ 



2n^ , , 

— y = 0. 



(A13) 



This is the differential equation of Legendre functions 
which in general only has solutions that are singular at 
^ = ±1 pjl, that is at a; = ±i?. The only exceptions are 



20? 



n{n -I- 1), 



(A14) 



where n is an integer. In this case the well-known Leg- 
endre polynomials Q 



(A15) 



solve Eq. (A13). Furthermore, they fulfill the orthogo- 
nality relation 



-1-1 

J PniOPrniOdC^O for n ^ TO. 



(A16) 



The frequencies of the elementary excitations are there- 
fore given by Eq. (|A14| ). The solution of the homoge- 
neous part of Eq. ( [A10| ) reads 

oo 

6p{x,t) ^"^AnSm{nnt + ipn)Pnix/R), (A17) 

71=1 

where An and tpn follow from the initial deviation from 
the stationary solution. 

Equation (A17) shows that the excitations do not de- 
cay, even in the presence of pump and loss terms in our 
generalized Gross-Pitaevskii equation, Eq. ([l7|). How- 
ever, these excitations do not grow but ra ther oscillate. 
Our numerical solution of Eq. (^7|) in Sec. [I C confirms 
this result. 

In Fig . p](b ) we do not find all frequencies allowed 
by Eq . ( |A14 ) since in our numerical solution of Eqs. 
( |l7| ) and (|18) we start with a symmetric initial condi- 
tion ^{x,t = 0). Because Eqs. ( p7| ) and ( [l8[ ) do not 
destroy this symmetry, p{x, t) as well as Sp{x^, the de- 
viation from the (symmetric) stationary solution of the 
generalized GPE, will always be symmetric. We therefore 
can only find excitations which correspond to Legendre 
polynomials of even order if we start with a symmetric 
tpix, t = 0). The corresponding frequencies are 



^2n ~ '^\/ n{2n + 1). 



(A18) 



This is also true for an anti-symmetric initial wave func- 
tion since the corresponding density is symmetric. 

Two other facts are worth mentioning: The excitations 
discussed above do not change the number of atoms. Us- 
ing Po(C) = 1 we find with the help of the orthogonality 
relation Eq. ( |A16D 



Sp{x, t) dx — 0. 



The second moment 



x^{t) 



^ x'^\ip{x,t)\'^ dx 
/|V(a;,t)|2dx 



can only oscillate with the frequency ^2 
follows from the relation 



e^lp2{o + \p^{o 



(A19) 

(A20) 
V^w. This 

(A21) 



together with the orthogonality relation, Eq. ( Al6| ). A 
generalization is the statement that 
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Pk{xlR) = 



J Pkix/R)\ijix,t)\^dx 



only shows oscillations with frequency Qk- 



(A22) 
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